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Abstract 

We propose a new criterion to analyse the order of phase transitions within 
a finite size scaling analysis. It refers to response functions like order pa- 
rameter susceptibilities and the specific heat and states different monotony 
behaviour in volume for first and second order transitions close to the transi- 
tion point. The criterion applies to analytical and numerical studies of phase 
diagrams including tricritical behaviour. 
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1 Introduction 



Finite size scaling analysis is particularly useful in numerical simulations to distin- 
guish first from second order phase transitions by their different onset in a large but 
finite volume. Usually it is the specific scaling behaviour of the peak in a suscep- 
tibility, the shift in the critical coupling/temperature and the width of the critical 
region that anticipates the nature of the transition in the thermodynamic limit [![. 
In this letter we formulate a new criterion in a finite volume to distinguish first from 
second order phase transitions. For a certain interval of the scaling region response 
functions with a nonanalytic behaviour in the infinite volume limit show different 
monotony behaviour for 1st and 2nd order transitions. Examples for such functions 
are the specific heat and order parameter susceptibilities. They are increasing in 
volume in a certain neighbourhood of T c for 2nd order transitions, and decreasing 
for 1st order transitions for some range in the scaling region, which will be specified 
below. 

Figure 1: (t, v )-plane in the vicinity of a phase transition at t — v — 0. t denotes the 
scaling field t = (T — T C (L)) /T C (L) , v is inverse to some power of the volume L x with 
x > 0. Let us consider the susceptibility as a function of t and v. Then, for a 1st order 
transition, x(t, v i) > x(t> v 2), whereas for a 2nd order transition x(t, v i) < x(ti v 2)- 
The shaded part is the crossover region between different asymptotic behaviour as 
(*,t;)-0. 




t 



For definiteness we fix the notation in terms of an order parameter susceptibility 
x(T, L), considered as a function of the temperature T and the spatial size parameter 
L, so that a symmetric D-dimensional volume is of size L D . By T C (L) we denote 
the location of the maximum of x(',L) in volume L D associated with the phase 
transition at T c = T c (oo). In the infinite volume 

x(t + T c (oo), L = oo) < oo as t — > (1) 

for a 1st order transition at T c , with a possible discontinuity at t = in the associated 
order parameter, t denoting the scaling field t — (T — T C (L))/T C (L), whereas 

X(t + T c (oo),L = oo) ~ ^|t|" 7 ast^O (2) 
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for a 2nd order transition with critical exponent 7 > 0. On the other hand, in the 
opposite order of limits, x(T c (L), L) diverges in both cases as L approaches infinity. 
Stated differently, at T C (L), in the limit L — > 00 x nas a 5-function singularity or 
a power law type of singularity for a 1st or 2nd order transition, respectively. It is 
this difference that is responsible for the different monotony properties for t 7^ in 
the finite volume. 

The essential statement of the monotony criterion is the following (details will 
be given below): For a sufficiently large size L s < 00, there is always a size L\ with 
L\ > L s and such that (cf. Fig. [I]) 

X (5 + T C (L), L > U) < X (5 + T C (L S ), L s ) for 1st order, (3) 
x{5 + T c (L),L>L l )>x{5 + T c (L s ) 1 L s ) for 2nd order, (4) 

that is, x(<5 + T C (L), L) is decreasing or increasing in volume for a 1st or 2nd order 
phase transition, respectively. The susceptibility is measured at fixed distance 5 to 
the volume dependent maximum, constraint by 

c x a(L)^ < \8\ < c 2 a(L s ). (5) 

Here cr(L) = L~ D or cr(L) = L^ 1 ^ for 1st or 2nd order transition, respectively, with 
v being the critical exponent of the correlation length, £ ~ |T — T c \~ u as T — > T c . 
Hence cr(L) is the width of the critical region in the volume L D , c\, c 2 and e are 
positive constants, typically e = 1. Beyond the general constraint that both linear 
sizes L s and L\ have to be sufficiently large, in addition L\ has to be sufficiently 
larger than L s , so that <y{L{) is considerably smaller than a(L s ). The monotony 
behaviour of x does not refer to values of Li close to L s . In terms of Fig. |TJ, L s and 
L\ are separated by the shaded area, for which we do not make any predictions on 
the behaviour of x- 

In section 2 we make these statements more precise. We specify the conditions 
under which the criterion holds and give the proof. As it turns out, the criterion 
applies whenever the standard conjectures on finite size scaling behaviour hold. 
Rigorous proofs of these standard conjectures are missing in general, neither do we 
attempt to give such proofs in the treatment below. In section 3 we conclude with 
some remarks on testing the criterion in analytical and numerical calculations. 

2 Statements of the criterion 

Let t denote the scaling field, i.e. t = (T—T C (L))/T C (L). T C (L) locates the maximum 
of the susceptibility in the volume L D . Furthermore we set v = L~ m with some 
m > 0. The infinite volume limit is obtained as v — > from above. 
The transition region is given by small t and v. Let 

H 2 := {(t,v) G R 2 I v > 0} (6) 

denote a half plane, U C H 2 the intersection of an open neighbourhood of G R 2 
with H 2 and W = U \ {0}. 
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2.1 First order transitions 



The typical behaviour of a susceptibility close to a first order transition is summa- 
rized by the following definition (with v = L~ D ). As usual, C{U*) denotes the set 
of real valued continuous functions in U* . 

Definition 1 Letu > 0. ^f(U) is the set of functions x G C(U*) with the following 
properties. 

1. x £ C X {U* \ {(R, 0)}), that is, x ^ continuously differentiable for v ^ 0. 

2. There exists u(t) > such that for all v < \t\/i>(t) 

I x(t,v) | < u. 

3. With appropriate positive constants c, Ki, K 2 and e we have in U* for v ^ 



V V 6 

I alX(*,w) | < 

at f z 

It follows already from this definition that v(t) is bounded from below by some 
positive constant. This is part of the following proposition. The main statement is 
the monotony behaviour of a function that obeys Definition |l|. 

Proposition 1 Let uj > 0, x ^(W). There exist positive 5, e such that the 
following statements hold. 

1. For all t with \t\ < e5, 

u(t) > e. (7) 

2. For all (t, v), (t, w) G U* with v{t)w < \t\ < ev and v < 5, we have w < v and 

X (t,v) > \ X (t,w)\. (8) 



The proposition includes the case w = of infinite volume. In this case, in the 
critical region \t\/v < e, the susceptibility x is always larger in the finite volume 
than in the infinite one. 

Proof: Let xe^(W),w>0. 

a. There are e, S > such that for all t, v with v < 8, \t\ < ev 

X (t,v)>2u. (9) 
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For, 

X(t,v) = x(0,w) 
and there are c , K > so that 



f 1 d 

7o ds d^ xivX 



r)=st 



x(0,iO>-, 

V 



K 

< — 

ir 



in U*. Hence, with e = c /(2K), 



*(*,«)>-- 



> c 



Co 

2u' 



V" V 

Choosing 6 = co/(4u;), we get for v < S the lower bound (|9]) on v). 
(3. By definition of v{t), we have whenever v{t)w < \t\ the upper bound 

\ X (t,w)\<cu. (10) 

7. Summarizing a and /?, it follows that for all w, v with u{t)w < \t\ < ev , v < 5, 
the strict bound 

\x(t, w )\ < X(t,v) 

holds. Furthermore, ( |10|) implies that \x{t,w')\ < u for all < w' < uo. Hence it 
must hold that w < v. It follows that 



t 



\t\ 



and hence v{t) > e for \t\ < e5. □ 

The remainder of this section is devoted to the function class As an 

example of such a function we consider the following typical representation of the 
magnetic susceptibility in the volume L D close to a first order phase transition, 



X2 (T, L) = c 2 L D exp (-f 2 L 2U (T - T c {L)f) + m (T, L). 



2D 1 



(11) 



Here c 2 , f 2 are positive real constants and t] 2 (-,L) is analytic for L < 00, locally 
uniformly convergent as L f 00 (so that f] 2 (-, 00) is analytic). Furthermore, T C (L) = 
T c - d 2 L~ m + o(L-' m ), with m > 0. 

We can view \ 2 as a function of some in essentially two ways. The first one 
corresponds to the point of view given in the introduction, in which the scaling field 
t measures the distance of the temperature T to the volume dependent location of 
the peak T C (L). It does not refer to a particular volume dependence of T C (L) itself. 
That is, m does not need further specification except that m > 0. 



We set v = L , t = (T - T C (L))/T C (L) and 



where r](t, v) 



X{t,v) = X2(T,L) 
= V2(T,L) and f(v) 



cv 1 exp (— f{y )v 2 t 2 ) +r](t,v), 



(12) 



f 2 T c {L)\ c = c 2 . 
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Lemma 1 X G \P^(ZY) for some uo > and sufficiently small U. v(t) of Def. [1] 
satisfies ^(t)^ 1 — > as |£| — >■ 0. 

Proof: For a, /3 > and a; G R the function 

/a,/? (a) = x " exp(-/3x 2 ) 
attains its maximum at x 2 = a/ (2/3), is positive and monotone decreasing for x > 

1 /2 

(a/ (2(3)) , and satisfies on R the bound 



a 
We 



We prove that t>) match the criteria of Definition |T| to belong to some ty±(U). 
There is 5 > such that for all v < 5 we have / < f(v) < / for some /, / > 0. We 
choose W such that (t,v) G W implies that v < 8 and |£| < f 1 / 2 and put 



^ = max (\r](t,v)\, \d t r](t,v)\) . 



Condition 1. is obviously satisfied. To show 2. we notice that for 



^ > v(t) > 
v \2f, 



one obtains 

\x(t,v)\ < 



c \t\ 

W\ ~v 



exp -f(v) 











+ 7] 


< c 









v(t) 



exp(-/(t;)z/(t) 2 ) + |r ? |. 



We choose u(t) = \t\ 1 and obtain with to — c/ (fe) + K v 

\x(t,v)\ < u. 
To verify the 3rd condition, we compute 

\ X {0,v)-cv- 1 \ = \n(t,v)\<K n , 
and notice that there is K > such that 



< 



2\t\ 

v 2 

V 



< -±rf\x\ + \d t v(t,v)\ 
l 



2fe 



+ f*K v 



+ K n 



K 

v 



Hence, x e *i(W). □ 

The next lemma states that the susceptibility %2 matches the criteria of Defini- 
tion ^ also if the scaling field t measures the distance of T from the phase transition 
point, i.e. to T C (L = oo). The only additional constraint then is that the peak T C (L) 
of X2 converges sufficiently fast to T c (oo), namely that 

\T C {L) -T c (oo)| < withm>D. (13) 

We set v = L~ D , t = (T - T c (oo))/T c (oo) and 

X (t,v) = X 2(T,L) =cv~ 1 exp(-fv- 2 (t + dv n ) 2 ) +7]{t,v), (14) 

where now rj(t,v) = r] 2 (T,L) and / = / 2 7c(oo) 2 , d = d 2 /T c (oo), and n = m/D. 
Similarly as above, we get 

Lemma 2 Let n > 1. TTten % G (W) /or some > and sufficiently small U. 
u(t) of Def. satisfies v{t)- 1 ^ as \t\ -> 0. 



Proof- For a, 6, e G R, e 7^ 0, 



We choose 5, 77 > sufficiently small such that the following estimates on x are true 
for all t, v with \t\ < 77, < v < 8, identifying U. 

X obviously satisfies condition 1. of Def. [l|. To show 2. we estimate 



\ X {t,v)\ < cv- 1 exp(-fv~ 2 (±t 2 -(dv n ) 2 )) + 

f 

< 2cv~ 1 exp (— — v~ 2 t 2 ) + \q\ . 



(15) 



Proceeding as in the proof of the previous lemma we obtain condition 2, with ap- 
propriate to. To realize the 3rd condition, we first notice that 



|x(0, v ) — cv I < cv 



-1 



+ \v\<— t 



for some positive constant K, with e = min (l,2(n — 1)) > 0. Finally there is K' > 
such that 



2fir 2 \t + dv n \\x\ + \d t r](t,v)\ 



< 2fcv- 



2/e 



K' 

• 2 = — • 



(16) 



Hence x e 



□ 
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Finally, we state without proof that every function x '■ — > R of the form 

Xtt,v) = -f(-) + x(t,v) (17) 

V V 

belongs to ^ (U) for some U and with appropriate ui > 0, if the following conditions 
are satisfied. 

la. x(t,v)eC 1 (W), 

lb. x(t, v) together with its (first) partial derivatives are uniformly bounded in U* . 

2a. / G C 1 (R) is a nonnegative function with /(0) > 0, 

2b. lim^^-i-oo \x\ l+e f(x) = for some e > 0, 

2c. (d/dx)f(x) is uniformly bounded on R. 

Any such function has the property that it "approaches 5" locally, i.e. there is 
c > so that for any test function g{t) supported in a small neighbourhood of t — 

lim i Mt, v)]g = lim / dt x(t,v)g(t) = cg(0). 
v— >o+ t)^o+ y 

The example discussed above corresponds to 

g(x) = (-) 1/2 exp (-ex 2 ) , c> 0. 

7T 

2.2 Second order transitions 

In contrast to the first order case, at a second order transition the susceptibility x can 
be divergent in the infinite volume limit as the critical temperature is approached. (It 
need not be divergent, since another second derivative of a suitable thermodynamic 
potential may be divergent instead.) If x is divergent with a critical exponent, it 
is described by the following definition (with v = L" 1 ^, v > being the critical 
exponent of the correlation length £). 

Definition 2 Let 7 > 0. ^(U) is the set of functions x £ C(U*) that satisfy the 
following conditions. 

1. There are constants A, K , e > such that in U* 

I x (t,0) - A |t|" 7 I < K \t\~ 1+e . 

Furthermore, there exist u, C > so that for \t\ > v v, 

X (t,v) > C X (t,0). 
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2. There are constants rj, B > such that for \t\ < r] v 

| X (t,v) | < Bv~\ 

Functions that obey Definition ^| are always increasing in volume close to the 
transition. 

Proposition 2 Let 7 > and x £ ^ZC^O- There are constants v , e > such that 
for all t, v, w with uw < \t\ < ev we have w < v and 

I X(t,v) I < x(t,w). 

The inequality includes the case w = 0, corresponding to infinite volume. That 
is, the susceptibility is always smaller in finite volume than in the infinite volume 
as long as we are in the critical region \t\/v < e. 

Proof: Let x £ ^K^Oj 7 > 0. There are numbers C,T>, v > such that in U* 
for |i| > v w 

X (t,w) > c x (t,o) > v \t\-\ 

Furthermore, there are r], B > such that for \t\ < r\ v 

\ X (t,v)\ < Bv-\ 

We choose e = min (r), {V /B) 1 ^ 1 ) and get for v w < \t\ < e v 

( e \ 7 Be 1 
\x(t,v)\ < Bv 7 < Bl — \ < — X (t,w) < x(t,w). 

This proves the lemma. □ 

The specific property for a 2nd order transition that the singular part of the 
free energy density behaves as a generalized homogeneous function implies for the 
susceptibility in a volume L D a typical form like 

X2 (T, L) = \T — T c r Q((T - T C )LV") + r](T, L), (18) 

with some 7 > 0. Here rj(-,L) has similar analyticity properties as in above, 
v > is the critical exponent of the correlation length £ ~ \T— T c \~ u . Q is continuous 
and behaves as 

lim |x|" 7 Q(x) = K > 0, 
lim Q(x) = C>0. 

x— >±oo 

The first equation expresses the absence of a divergence of xi f° r finite L, the second 
one its presence in the infinite volume case. With t = (T — T c )/T c , v = L~ l l u and 

X (t,v) := X 2(T,L) 
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it is straightforward to verify that x belongs to ^1(U). 
More generally, every function x : W* — > R of the form 

X(t,v) = 1/(1) +*(*,„) (19) 

with 7 > belongs to some ^(U), if the following conditions are satisfied: 
1. x(t,v)eC 1 (W). 

2a. / G C\B), and /(0) > 0, 

2b. lim^.^-1-oo |x| 7 /(x) = C for some finite C > 0. 

Compared to the first order case fll7|) the essential difference comes from property 
(2b). 

As an example, 

X (t,v) = {f + v 2 )- (m/2) ,m>0, 

belongs to the class 

3 Applications 

The monotony criterion has been applied in the framework of convergent series 
expansions of O(N) symmetric spin models on the lattice, performed to the 20th 
order in the hopping parameter (corresponding to the inverse temperature) both in 
the infinite and finite volume 0. By means of the series representation it is hard 
to compute the scaling behaviour of a susceptibility x(T, L) at the maximum T C (L), 
since T C (L) as well as the transition at T c (oo) are extrapolated as the convergence 
radius of the series. The monotony criterion circumvents this difficulty in that it 
allows us to determine the order of the transition by the different L-dependence of 
x{T,L) already close to T C (L), that is by an increase in volume for a second order 
transition, and a decrease for a first order transition. 

In series representations the larger volume Lf is conveniently chosen to be in- 
finite so that df)-(||), with L = oo and cr(L) = 0, give strong criteria in the whole 
scaling region. In Monte Carlo simulations both volumes Lf and Lf must be finite. 
The monotony criterion works equally well in this case. Both volumes have to be 
sufficiently large in order to match the conditions of the criterion (which reflect the 
standard assumptions of finite size scaling). Some care is needed to ensure that the 
susceptibility is actually measured at a temperature that satisfies (||p, i.e. out of 
a region where the difference in the volume dependence between first and second 
order transitions appears. In particular this concerns the lower bound on 5 in (^) 
because for L\ < oo, also in case of a first order transition, a small neighbourhood 
around the peak at T C (L) exists, where x(T, L) is increasing in volume because of 
the rounding of the 5-singularity. 
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To summarize, let A denote a generic coupling constant. Application of the 
monotony criterion amounts to a choice of two sizes L s and Li with L s , Li and 
Li/L a sufficiently large. Let us define the ratio r 



r(L s ,Li;X) := 1 



X(T C (L S ) + 5,L S ;X) 
X (T C (L0 + 5,L Z ;A)' 



(20) 



with 5 satisfying the constraint ([5]). The monotony criterion then says that 



In (|20|) the susceptibility is measured at fixed distance to the volume-dependent 
location of its maximum. A reasonable alternative is to measure x{T, L) for both 
volumes at the same temperature T regardless of L, but T chosen as T = T c (oo) + 5. 
This requires that T C (L) approaches T c = T c (oo) according to 



cf. @ and (|T3|). The behaviour ( pT|) meets the standard finite size scaling behaviour. 
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> , 2nd order 
< , 1st order 

= , tricritical point for r changing sign at A. 



T c (L)-T c (oc)| < ca(L) 



(21) 
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